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Q ■ Abstract 

CO ' In this paper, the concept of sparse difference resultant for a transformally essential 

system of difference polynomials is introduced and its properties are proved. In partic- 
ular, order and degree bounds for sparse difference resultant are given. Based on these 
. bounds, an algorithm to compute the sparse difference resultant is proposed, which is 

' single exponential in terms of the number of variables, the Jacobi number, and the size 

, of the transformally essential system. Also, the precise order, degree, a determinant 

O ' representation, and a Poisson-type product formula for difference resultants are given. 

Keywords. Sparse difference resultant, difference resultant, Laurent transformally es- 
^ ' sential system, Jacobi number, single exponential algorithm. 
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. 1 Introduction 

cn 

CN . The resultant, which gives conditions for an over-determined system of polynomial equations 

to have common solutions, is a basic concept in algebraic geometry and a powerful tool in 
elimination theory O (Tj [U dTJ [l8l [26] . The concept of sparse resultant is originated from the 
work of Gelfand, Kapranov, and Zelevinsky on generalized hyper geometric functions, where 
the central concept of ^-discriminant is studied [16]. Kapranov, Sturmfels, and Zelevinsky 
introduced the concept of ^-resultant [19j. Sturmfels further introduced the general mixed 
■ sparse resultant and gave a single exponential algorithm to compute the sparse resultant 

|26t |27] . Canny and Emiris showed that the sparse resultant is a factor of the determinant of 
a Macaulay style matrix and gave an efficient algorithm to compute the sparse resultant based 
on this matrix representation [Uj . A determinant representation for the sparse resultant was 
given by D 'Andrea [8]. Recently, in [T3], a rigorous definition for the differential resultant of 
n + 1 generic differential polynomials in n variables was presented [H] and also the theory 
of sparse differential resultants for Laurent differentially essential systems was developed 
\22\ |23] . It is meaningful to generalize the theory of sparse resultant to difference polynomial 
systems. 

In this paper, the concept of sparse difference resultant for a Laurent transformally 
essential system consisting of n -|- 1 Laurent difference polynomials in n difference variables 
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is introduced and its basic properties are proved. In particular, we give order and degree 
bounds for the sparse difference resultant. Based on these bounds, we give an algorithm 
to compute the sparse difference resultant. The complexity of the algorithm in the worst 
case is single exponential of the form 0(m'^("''"^ ^(nJ)'^*^'"^)), where n,m,J, and I are the 
number of variables, the degree, the Jacobi number, and the size of the Laurent transformally 
essential system respectively. Besides these, the difference resultant is introduced and its 
basic properties are given, such as its precise order, degree, determinant representation, and 
Poisson-type product formula. 

Although most properties for sparse difference resultants and difference resultants are 
similar to its differential counterpart given in [22\ [23l [T3] , some of them are quite different 
in terms of descriptions and proofs. Firstly, the definition for difference resultant is more 
subtle than the differential case as illustrated by Problem 13.151 in this paper. Secondly, the 
criterion for transformally essential systems given in Section 3.3 is quite different and much 
simpler than its differential counterpart given in [23]. Also, a determinant representation for 
the difference resultant is given in Section [6l but such a representation is still not known for 
differential resultants [30^ I25j. Finally, some properties are more difficult in the difference 
case. For instance, we can only show that the vanishing of the difference resultant is a 
necessary condition for the corresponding difference polynomial system to have a common 
nonzero solution. While, the sufficient condition part is still open. Also, there does not exist 
a definition for homogeneous difference polynomials, and the definition we give in this paper 
is different from its differential counterpart. 

The rest of the paper is organized as follows. In Section O we prove some preliminary 
results. In Section [3l we first introduce the concepts of Laurent difference polynomials and 
Laurent transformally essential systems, and then define the sparse difference resultant for 
Laurent transformally essential systems. Then basic properties of sparse difference resultant 
are proved in Section HI And in Section [5l we present an algorithm to compute the sparse 
difference resultant. Then we introduce the notion of difference resultant and give its basic 
properties in section [6l In Section [71 we conclude the paper by proposing several problems 
for future research. 

2 Preliminaries 

In this section, some basic notations and preliminary results in difference algebra will be 
given. For more details about difference algebra, please refer to [5l [21] . 

2.1 Difference polynomial ring 

An ordinary difference field J-" is a field with a third unitary operation a satisfying that 
for any a,b £ T, a{a + b) = a{a) + cr(6), a{ab) = a{a)a{b) and a{a) = if and only if 
a = 0. We call a the transforming operator of If a G J^, (T{a) is called the transform of 
a and is denoted by a^^\ And for n € Z"*", a"'{a) = a'^~^{a{a)) is called the n-th transform 
of a and denoted by a^"'\ with the usual assumption By a'"] we mean the set 

{a, }. A typical example of difference field is Q(x) with a{f{x)) = f{x + 1). 

Let S" be a subset of a difference field G which contains J-. We will denote respectively 
by T[S], J^{S), J^{S}, and T{S) the smallest subring, the smallest subfield, the smallest 
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difference subring, and the smallest difference subfield of G containing T and S. If we 
denote G(5) = {a''a\k > 0,a € S}, then we have J^{S} = T[Q{S)] and T{S) = T{e{S)). 

A subset 5 of a difference extension field ^ of is said to be transformally dependent 
over if the set {(7*^010 G 5,/c > 0} is algebraically dependent over and is said to be 
transformally independent over J^, or to be a family of difference indeterminates over in 
the contrary case. In the case S consists of one element a, we say that a is transformally 
algebraic or transformally transcendental over T respectively. The maximal subset il. G 
which are transformally independent over J-" is said to be a transformal transcendence basis 
of Q over J^. We use Atr.deg^/J^ to denote the difference transcendence degree of G over 

which is the cardinal number of $7. Considering T and G as ordinary algebraic fields, we 
denote the algebraic transcendence degree of G over T by ti.degG / J^- 

Now suppose Y = {yi,y2, ■ ■ ■ ,yn} is a set of difference indeterminates over T. The 
elements of J^{Y} = J-[y^p : j = l,...,n;k E Nq] are called difference polynomials over 
J-" in Y, and J-'{Y} itself is called the difference polynomial ring over J-" in Y. A difference 
polynomial ideal I in J^{Y} is an ordinary algebraic ideal which is closed under transforming, 
i.e. <t(I) C I. If I also has the property that a^^^ G I implies that a £ I, it is called a 
reflexive difference ideal. And a prime (resp. radical) difference ideal is a difference ideal 
which is prime (resp. radical) as an ordinary algebraic polynomial ideal. For convenience, a 
prime difference ideal is assumed not to be the unit ideal in this paper. If S is a finite set of 
difference polynomials, we use (S) and [S] to denote the algebraic ideal and the difference 
ideal in /"{Y} generated by S. 

An n-tuple over T is an n-tuple of the form a = (oi, . . . , a„) where the Oj are selected 
from some difference overfield of J-. For a difference polynomial / G . . . , a is 

called a difference zero of / if when substituting y^ by a- in /, the result is 0. An n-tuple 
T] is called a generic zero of a difference ideal I C J^{Y} if for any polynomial P G J-'{Y} we 
have P(ry) = <^=> P G X. It is well known that 

Lemma 2.1 ^-77/ A difference ideal possesses a generic zero if and only if it is a reflexive 
prime difference ideal other than the unit ideal. 

Let X be a reflexive prime difference ideal and ij a generic point of Z. The dimension of 
Z is defined to be Atr.degJ^(77)/J^. 

Given two n-tuples a = (oi, . . . , a„) and a = (ai, . . . , a„) over T. a is called a specializa- 
tion of a over X", or a specializes to a, if for any difference polynomial P G J-'{Y}, -P(a) = 
implies that P{sl) = 0. The following property about difference specialization will be needed 
in this paper. 

Lemma 2.2 Let Pi(U,Y) G J'(Y){U} {i = l,...,m) where U = {ui,...,Ur) and Y = 
(yi, . . . ,y„) are sets of difference indeterminates. //Pj(U, Y) {i = 1, . . . ,m) are transformally 
dependent over J^{V) , then for any difference specialization U to U which are elements in T , 
Pi(U, Y) (z = 1, . . . , m) are transformally dependent over T . 

Proof: It suffices to show the case r = 1. Denote u = u\. Since Pi{u,Y) {i = l,...,m) 
are transformally dependent over J-{u), there exist natural numbers s and I such that 
Ff'\u,Y) {k < s) are algebraically dependent over T{u^^^\k < s + I). When u special- 
izes to u & u^'^^ {k > 0) are correspondingly algebraically specialized to u^^^ G J-". By [291 
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p. 161], P- (n, Y) [k < s) are algebraically dependent over J^. Thus, Pi{u,Y) (i = 1, . . . ,m) 
are transformally dependent over T. □ 



2.2 Characteristic sets for a difference polynomial system 

Let / be a difference polynomial in J-'{Y}. The order of / w.r.t. yi is defined to be the 
greatest number k such that y^''^ appears effectively in /, denoted by ord(/, yj). And if 
yi does not appear in /, then we set ord(/, yj) = —oo. The order of / is defined to be 
maxjord(/,yi), that is, ord(/) = max, ord(/, yi). 

A ranking ^ is a total order over G(Y) = {a^yi\l < i < n,k > 0}, which satisfies the 
following properties: 

1) a{e) > e for all derivatives 9 G e(Y). 

2) 9i>02^ cj(^i) > a{92) for 0i,02 G e(Y). 

Let / be a difference polynomial in J^{Y} and ^ a ranking endowed on it. The greatest 
yj^^ w.r.t. ^ which appears effectively in / is called the leader of p, denoted by ld(/) and 
correspondingly yj is called the leading variable of /, denoted by lvar(/) = yj. Let the degree 
of / in ld(/) be d. The leading coefficient of / as a univariate polynomial in ld(/) is called 
the initial of / and is denoted by Ij . 

Let p and q be two difference polynomials in J^{Y}. q is said to be of higher rank than 

p if 

1) ld{q) > ld(p), or 

2) ld{q) = ld(p) = yf^ and deg{q,yf^) > deg{p,yf^). 

Suppose ld{p) = ■ q is said to be reduced w.r.t. p if deg((7, yj'^^''') < deg(p, yj'^^) for 
all I G No. 

A finite chain of nonzero difference polynomials A = Ai, . . . , Am is said to be an ascending 
chain if 

1) m = 1 and yli 7^ or 

2) m > 1, Aj > Ai and Aj is reduced w.r.t. Ai for 1 < i < j < m. 

Let A = Ai, A2, ■ ■ ■ , At be an ascending chain with Ij as the initial of Ai, and / any 
difference polynomial. Then there exists an algorithm, which reduces / w.r.t. ^ to a 
polynomial r that is reduced w.r.t. A, satisfying the relation 

flJ](a^I,r--/^r,mod [A], 

i=l k=0 

where the Cj^ are nonnegative integers. The difference polynomial r is called the difference 
remainder of / w.r.t. A |15j . 

Let A be an ascending chain. Denote to be the minimal multiplicative set containing 
the initials of elements of A and their transforms. The saturation ideal of A is defined to be 

sat(^) = [A] : Ia = {p ■■ e Ia, s.t. hp £ [A]}. 

And the algebraic saturation ideal of A is asat(^) = (A) : I^, where is the minimal 
multiplicative set containing the initials of elements of A. 
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An ascending chain C contained in a difference polynomial set S is said to be a character- 
istic set of S, if S does not contain any nonzero element reduced w.r.t. C. A characteristic 
set C of a difference ideal J' reduces to zero all elements of J'. 

Let ^ be a characteristic set of a reflexive prime difference ideal X. We rewrite A as the 
following form 

{^11, • • • ) ^ifci 
... 
Api , . . . , ^pfcp 

where lvar(j4jj) = 1/^ for j = 1, . . . , ki and ord(Ajj, ycj < ord(Aji, ycj for j < I. In terms of 
the characteristic set of the above form, p is equal to the codimension oil, that is n — dim(X). 
Unlike the differential case, here even though I is of codimension one, there may be more 
than one difference polynomials in a characteristic set of Z as shown by the following example. 



Example 2.3 Let An = + + 1, A12 = yf^ — yi. ThenI = [An, ^12] is a reflexive 

prime difference ideal whose characteristic set is A = Aii,Ai2 and I = sat(^) fT5^ . Note 

(2) (2) 

that [All] is not a prime difference ideal, because a{Aii) — An = {y\ —yi)(y\ +yi) € [^11] 

(2) (2) 
and both yl — yi and yl + yi are not in [An]. 

Now we proceed to show that a property of uniqueness still exists in characteristic sets 
of a reflexive prime difference ideal in some sense. Firstly, we need several algebraic results. 

Let B = Bi, . . . , Bm be an algebraic triangular set in T[xi, . . . , Xn] with lvar(i?j) = yt 
and U = {xi, . . . , Xn}\{yi, ■ ■ ■ ,ym}- A polynomial / is said to be invertible w.r.t. A if 
{f,Ai,..., As) n K[U] 7^ {0} where lvar(/) = Ivar(As). We call B a regular chain if for each 
i > 1, the initial of Bi is invertible w.r.t. Bi, . . . , -Bi_i. For a regular chain B, we say that / 
is invertible w.r.t. asat(^) if {f, asat{B)) n T[U] ^ {0}. 

Lemma 2.4 Let B be a regular chain in J-[xi, . . . , Xn]. If \l asat(;B) = Hi^i "^j irre- 
dundant prime decomposition of y/asatJW), then a polynomial f is invertible w.r.t. asat(;S) 
if and only if f ^ Vi for all i = 1, . . . ,m. 

Proof: Since asat(S) = fli^i is an irredundant prime decomposition of asat(S), U 
is a parametric set of Vi for each i by [13j . And for prime ideals Vi, f ^ Vi \i and only 
if {f,Vi) r\J'[U] / {0}. If / is invertible w.r.t. asat(6), {0} / (/,asat(S)) n /"[^J] C 
1^ ^\U]- Thus, f ^ Vi for each i. For the other side, suppose f ^ Vi for all i, then 
there exist nonzero polynomials hi{U) such that hi{U) S {f,Vi). Thus, there exists t gN 
such that (Hi^i hi{U)Y G (/, asat(0)). So / is invertible w.r.t. asat(;B). □ 

Lemma 2.5 J^l Let B be a regular chain in J-[U,Y]. Let f be a polynomial in T[U,Y] and 
L in J='[U]\{0} such that Lf E (B). Then f € asat(^). 

Lemma 2.6 Let A be an irreducible difference polynomial in J^{Y} with deg{A,yig) > for 
some io. If f is invertible w.r. t. = A,A(i),...,A('=) under some ranking then cr{f) is 
invertible w.r.t. aI'^+^I = A, . . . , A^^"*"^-* . In particular, A^'''^ is a regular chain for any k >0. 
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Proof: Since as a difference ascending chain, A is coherent and proper irreducible, by Theo- 
rem 4.1 in [15], A is difference regular. As a consequence, ^['^^ is regular for any k >0. □ 

The following fact is needed to define sparse difference resultant. 

Lemma 2.7 Let I be a reflexive prime difference ideal of codimension one in J^{Y}. The 
first element in any characteristic set of X w.r.t. any ranking, when taken irreducible, is 
unique up to a factor in T . 

Proof: Let A = A\^. . . , be a characteristic set of T w.r.t. some ranking 1% with A\ 
irreducible. Suppose Ivar(^) = y\. Given another characteristic set B = Bi,...,Bi of X 
w.r.t. some other ranking (Bi is irreducible), we need to show that there exists c (z J- 
such that Bi = c ■ Ai. It suffices to consider the case lvar(;B) 7^ yi. Suppose lvar(i?i) = 1/2 • 
Clearly, y2 appears effectively in Ai for B reduces Ai to 0. And since I is reflexive, there 
exists some zq such that deg{Ai,yi^^) > 0. 

Suppose ord(yli,y2) = 02- Take another ranking under which y^^^ is the leader of Ai 
and we use Ai to distinguish it from the Ai under By Lemma 12.61 for each k, A[ and 
are regular chains. 

Now we claim that asat{A^^^) = asat{A^^^). On the one hand, for any polynomial 
/ G asat(Afl), we have {Ui=o (^Ai)T f e (^f'). Since I^i is invertible w.r.t. Ai, by 
Lemma 12.61 is invertible w.r.t. aJ*'. Thus, (11^=0 ^H^Ai))" is invertible w.r.t. A^^\ 

Denote the parameters of A^^'^ by U. So there exists a nonzero polynomial h{IJ) such 
that h{U) G ((nLo^'(I^i))"'4''^)- Thus, h{U)f G (if^). Since if^ is a regular chain, 
by Lemma 12.51 / ^ asat{A^^^). So asat{Af^^) C asat(A^'^'). Similarly, we can show that 
asat(if') C asat(Af^). Thus, asat(Af^) = asat(if^). 

Suppose oid{Bi,y2) = o^- It is clear that 02 > o^- We now proceed to show that it 
is impossible for 02 > o^- Suppose the contrary, i.e. 02 > o'^- Then B\ is invertible w.r.t. 

asat(if^). Suppose ^asat(if') = flLi ^» 

is an irredundant prime decomposition. By 
Lemma 12.41 ^\ ^ T^i for each i. Since asat(Afl) = asat(if^), using Lemma 12.41 again, B\ 
is invertible w.r.t. a&a\,{^J^^\ Thus, there exists a nonzero difference polynomial B with 
ord(ff, yi) < ord(Ai,yi) such that B G (Si, asat(A5'^^)) C X, which is a contradiction. Thus, 
02 = 02. Since B reduces A\ to zero and A\ is irreducible, there exists c ^ T such that 
Bx = c-Ai. □ 



3 Sparse difference resultant 

In this section, the concepts of Laurent difference polynomials and transformally essential 
systems are first introduced, and then the sparse difference resultant for transformally es- 
sential systems is defined. And we also give a criterion for Laurent transformally essential 
systems in terms of the support of the given system. 
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3.1 Laurent difference polynomial 



Let J- be an ordinary difference field with a transforming operator a and J-'{Y} the ring 
of difference polynomials in the difference indeterminates Y = {yi, ...,?/„}. Similar to |23j . 
before defining sparse difference resultant, we first introduce the concept of Laurent difference 
polynomials. 

Definition 3.1 ^ Laurent difference monomial of order s is a Laurent monomial in vari- 
ables YH = (yf'^)i<j<„;o<fc<s- More precisely, it has the form nr=i 11^=0 (^i'^V*''' where dik 
are integers which can he negative. A Laurent difference polynomial over F is a finite linear 
combination of Laurent difference monomials with coefficients in T . 

Clearly, the collections of all Laurent difference polynomials form a commutative differ- 
ence ring under the obvious sum, product operations and the usual transforming operator 
CT, where all Laurent difference monomials are invertible. We denote the difference ring of 
Laurent difference polynomials with coefficients in T by T{y\^y'f^ ^ . . . ,yn,yn^}, or simply 



Definition 3.2 For every Laurent difference polynomial F € J-'{Y, Y^^}, there exists a 
unique laurent difference monomial M such that 1) M ■ F £ J^{Y} and 2) for any Laurent 
difference monomial T with T ■ F € ^{Y}, T ■ F is divisible by M ■ F as polynomials. This 
M ■ F is defined to he the norm form of F , denoted by ^{F). The order and degree of^{F) 
is defined to be the order and degree of F , denoted by ord(F) and deg(-F) . 

In the following, we consider zeros for Laurent difference polynomials. 

Definition 3.3 Let F he a Laurent difference polynomial in F{Y ,Y^^^ . An n-tuple (ai, 
. . . , On) over F is called a nonzero difference zero of F if for all i, ai and F{ai, . . . , a„) = 



For an ideal I € -FjYjY^^}, the difference zero set of I is the set of common nonzero 
difference zeros of all Laurent difference polynomials in I. We will see later in Example 14.2^ 
how nonzero difference solutions are naturally related with the sparse difference resultant. 

3.2 Definition of sparse difference resultant 

In this section, the definition of the sparse difference resultant will be given. Similar to 
the study of sparse differential resultants, we first define sparse difference resultants for 
Laurent difference polynomials whose coefficients are difference indeterminates. Then the 
sparse difference resultant for a given Laurent difference polynomial system with concrete 
coefficients is the value which the resultant in the generic case assumes for the given case. 

Suppose Ai = {Mio, Mil, . . . , Mil-} {i = 0,1,..., n) are finite sets of Laurent differ- 
ence monomials in Y. Consider n + 1 generic Laurent difference polynomials defined over 



by.F{Y,Y-i}. 



0. 



^0 ) • • • ! An ■ 
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where all the Uik are transformally independent over the rational number field Q. Denote 

Ui = {uio, Uii,..., Uin) {i = 0,...,n) and u = Ui^QUi\{uio} . (2) 

The number + 1 is called the size of Pj. To avoid the triviality, > 1 (i = 0, . . . , n) are 
always assumed in this paper. 

Definition 3.4 A set of Laurent difference polynomials of form ([Ip is called Laurent trans- 
formally essential if there exist ki (i = 0, . . . ,n) with 1 < ki < li such that Atr.degQ( , 

■ ■ ■ , ~^^)/Q. = ^- IiT- this case, we also say that Aq, • • • , An form a Laurent transfor- 
mally essential system. 

Although Mjo are used as denominators to define transformally essential system, the 
following lemma shows that the definition does not depend on the choices of Mjo- 

Lemma 3.5 The following two conditions are equivalent. 

1. There exist ki {i = 0, . . . ,n) with 1 < ki < li such that Atr.degQ(^^^, . . . , 



n. 



2. There exist pairs {ki,ji) (i = 0, . . . , n) with ki / € {0, . . . , li} such that Atr.deg ( 



'Mo 



'JO 



n. 



Proof: Similar to the proof of \23\ Lemma 3.7], it can be easily shown. □ 

Let m be the set of all difference monomials in Y and [N(Po), . . . ,N(P„)] the difference 
ideal generated by N(Pi) in Q{Y, uq, . . . , u„}. Let 

Iy,u = ([N(Po),...,N(P„)] :m). (3) 

The following result is a foundation for defining sparse difference resultants. 

Theorem 3.6 Let Po,...,Pn be Laurent difference polynomials defined in (Op. Then the 
following assertions hold. 

1. Xy,u is a reflexive prime difference ideal in Q{Y, uq, . . . , u„}. 

2. Zy,u n Q{uo, . . . ,u„} is of codimension one if and only if¥Q, . . . ,P„ form a Laurent 
transformally essential system. 

Proof: Let rj = (r/i, . . . ,r/„) be a sequence of transformally independent elements over Q(u), 
where u is defined in ([2]). Let 

^ Mik{ri) /. „ 1 ^ fA\ 



k=l 



We claim that 9 = (ry; Co, ^oi, . . . , uoio; ■ ■ ■ ; Cn,Uni, ■ ■ ■ , Uni„) is a generic point of Zy,u, which 
follows that Xy,u is a reflexive prime difference ideal. 
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Denote N(Pj) = MjPj (z = 0, . . . , n) where Mj are Laurent difference monomials. Clearly, 
N(Pj) = MjPj vanishes at 0(i = 0, . . . ,n). For any / G 2y,u) there exists an M € m such 
that Mf G [N(Po), . . . ,N(P„)]. It follows that f{9) = 0. Conversely, let / be any difference 
polynomial in Q{Y, uq, . . . , u„} satisfying f{e) = 0. Clearly, N(Po), N(Pi), . . . , N(P„) con- 
stitute an ascending chain with Uio as leaders. Let /i be the difference remainder of / w.r.t. 
this ascending chain. Then /i is free from Ujo {i = 0, . . . , n) and there exist a, s G N such 
that {U7=oUi=oi^\MiMio))r ■ f = /i, mod [N(Po), . . . , N(P„)]. Clearly, = 0. Since 

/i G Q{u,Y}, /i = 0. Thus, / G Xy,u. So Iy, u is a reflexive prime difference ideal with a 
generic point 6. 

Consequently, ly.u H Q{uo, . . . ,u„} is a reflexive prime difference ideal with a generic 
point C = {Co,uoi,- ■ ■ , uqi^;- ■ ■ ; Cn,Uni, ■ ■ ■ ,UniJ- FVom (HD, it is clear that Atr.degQ(C)/Q < 
^27=0 + If th^'^^ ^^ist P^i'^s (ikJk) (fc = 1, . . . , n) with I < jk < k^, and ik^ 7^ ik^ {ki / 
^2) such that -jT[^, • • • , ^"^g are transformally independent over Q, then by Lemma 



Cn!---)Ci„ are transformally independent over Q(u). It follows that Atr.deg( 
X]r=o + Thus, Iy,u n Q{uo, . . . , u^} is of codimension 1. 

Conversely, let us assume that Iy,u H Q{uo,...,u„} is of codimension 1. That is, 
Atr.degQ(C)/Q = '^^=oh + n'. We want to show that there exist pairs {ik,jk) {k = 1, . . . ,n) 
with 1 < jfc < and ik-^ 7^ ifcj {ki 7^ /C2) such that -jj^, ■ ■ ■ , are transformally inde- 



pendent over Q. Suppose the contrary, i.e., -^j^^^^, ■ ■ ■ , ^m"^^^^) transformally dependent 



transformally de] 



for any n different and jk G {1, . . . , / jj.}. Since each is a linear combination of ^ ^^^^ 

(ifc = li • • • 1 4fe)) it follows that Cn, • • • , Cin are transformally dependent over Q(u). Thus, we 
have Atr.degQ(C)/Q < Y17=o h + n, a contradiction to the hypothesis. □ 

Let [Pq, . . . , Pn] be the difference ideal in Q{Y, Y^-*^; uq, . . . , u„} generated by Pj. Then 
we have 

Corollary 3.7 [Pq, Pi, ... , P„] n Q{uo, . . . , u„} is a reflexive prime difference ideal of codi- 
mension one if and only if {Pj : i = 0, . . . ,n} is a Laurent transformally essential system. 

Proof It is easy to show that [Pq, Pi, . . . , P„] n Q{uo, . . . , u„} = Zy,u n Q{uo, . . . , u„}. And 
the result is a direct consequence of Theorem 13.61 □ 

Now suppose {Pq, . . . ,Pn} is a Laurent transformally essential system. Denote the dif- 
ference ideal [Pq, . . . , Pn] n Q{uo, . . . , u„} by I^- Since Xu is a reflexive prime difference ideal 
of codimension one, by Lemma 12.71 there exists a unique irreducible difference polynomial 
R(u; uoo, • • • , Uno) = R.(uo, . . . , u„) G Q{uo, . . . , u„} such that R can serve as the first poly- 
nomial in each characteristic set of 1^ w.r.t. any ranking endowed on uq, . . . , u„. That is, if 
Uio appears in R, then among all the difference polynomials in Xu, R is of minimal order in 
Uio and of minimal degree with the same order. 

Now the definition of sparse difference resultant is given as follows: 

Definition 3.8 The above R(uo, . . . , u„) G Q{uo, . . . , u„} is defined to be the sparse differ- 
ence resultant of the Laurent transformally essential system Pq, . . . , Pn, denoted by Res^^o^...^^^^ 
or RespQ^...^p„ . And when all the Ai are equal to the same A, we simply denote it by Res^^. 
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The following lemma gives another description of sparse difference resultant from the 
perspective of generic point, 

Lemma 3.9 Let Ci = ~ Ylk=i ^ifc M^o(^j — 0,1, ... ,n) defined as in equation Q), where 
r] = {rji, . . . ,rjn) is a generic point of [0] over Q(u). Then among all the polynomials in 
Q{uo, . . . , u„} vanishing at (u; Co, • • • , Cn), R-(uo, . . . , u„) = R(u; uoo, • • • , u„o) is of minimal 
order and degree in each Uio {i = 0, . . . ,n). 

Proof: It is a direct consequence of Theorem 13.61 and Definition 13.81 □ 



Remark 3.10 From its definition, the sparse difference resultant can be computed as follows. 
With the characteristic set method given in 115], we can compute a proper irreducible ascend- 
ing chain A which is a characteristic set for the difference polynomial system {¥q,¥i, . . . ,Pn} 
under a ranking such that Uij < y^- Then the first difference polynomial in A is the sparse 
difference resultant. This algorithm does not have a complexity analysis. In Section 5, we 
will give a single exponential algorithm to compute the sparse difference resultant. 

We give several examples to show sparse difference resultant. 



Example 3.11 Let n = 2 and Pj has the form 



(2) , (3) , (3) . . 



0,1,2) 



It is easy to show that yPVyi^^ '^^'^ V^'^ lUi''^ ^"^^ transformally independent over Q. Thus, 
Po,Pi,P2 form a Laurent transformally essential system. The sparse difference resultant is 



R 



^iOO ^^01 '"02 

uio nil ui2 

U20 U21 U22 



The following example shows that for a Laurent transformally essential system, its sparse 
difference resultant may not involve the coefficients of some Pj. 



Example 3.12 Let n = 2 and Pj has the form 

Po = lioo + uoiyiy2, Pi = uio + uiiy^^^y!^^\ P2 = U20 + ^212/2- 

Clearly, Po,Pi,P2 form a Laurent transformally essential system. And the sparse difference 
resultant o/Po,Pi,P2 is 

T3 (1) (1) 

R = ^00^11 - ^^01 ^10; 
which is free from the coefficients of P2 . 

The above example can be used to illustrate the difference between the differential 
and difference cases. If Po,Pi,P2 in Example 13.121 are differential polynomials, then the 
sparse differential resultant is UqiUioU2qU2i — uoiUqqUiiU2oU2iU2q + uooUqiUuU2oU2iU2o + 
noinoo'"iin2o(n2i)^ + uoo'Woiniin2i(n2o)^ - 2-uoi'Uoo'Uii 7x20^21^20^21 + uoiUqqUuuIqU2iU2i - 
uooUq-^^uiiU2iU2iU2q which contains coefficients of Po,Pi,P2. 
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Remark 3.13 When all the Ai{i = 0, . . . ,n) are sets of difference monomials, unless ex- 
plicitly mentioned, we always consider Pj as Laurent difference polynomials. But when we 
regard Pj as difference polynomials, Res^Q^...^^^ is also called the sparse difference resultant 
of the difference polynomials Pj and we call Pj a transformally essential system. In this 
paper, sometimes we regard Pj as difference polynomials where we will highlight it. 

We now define the sparse difference resultant for any set of specific Laurent difference 
polynomials over a Laurent transformally essential system. For any finite set A of Laurent 
difference monomials, we use C{A) to denote the set of all Laurent difference polynomials of 
the form ^a,/£_4«m-^ where the aj\,f are in some difference extension field of Q. Then C{A) 
can be considered as the set of all /-tuples over Q where / = \ A\. 

Definition 3.14 Let Ai = {Mio, Mn, . . . , Mu.} {i = 0,1,..., n) be a Laurent transfor- 
mally essential system. Consider n + 1 Laurent difference polynomials (Fq, Fi, . . . , F„) G 
nr=o'^("^«)' -^^^ sparse difference resultant of Fq, Fi, . . . , Fn, denoted as ResFo,...,F„) is ob- 
tained by replacing Ui by the corresponding coefficient vector of Fi in ^es_Ao,...,A„i'^Oj • • • ) 

A major unsolved problem about difference resultant is whether R defined above contains 
all the information about the elimination ideal 2y,u defined in ([3]). More precisely, we propose 
the following problem. 

Problem 3.15 As shown by Example \2.3[ the characteristic set for a reflexive prime differ- 
ence ideal could contain more than one elements. Let Iy,u be the ideal defined in (0). Then 
Iy,\i is a reflexive prime difference ideal of codimension one and 



where R is the sparse difference resultant of¥Q,...,¥n and R, . . . , Rm is a characteristic 
set o/Xy,u- We conjecture that m = 0, or equivalently Iy,u = sat(R), which is similar the 
differential case. If this is valid, then better properties can be shown for sparse difference 
resultant as we will explain later. 

3.3 Criterion for Laurent transformally essential systems in terms of the 



Let Ai{i = 0, . . . , n) be finite sets of Laurent difference monomials. According to Defini- 
tion 13. 4( in order to check whether they form a Laurent transformally essential system, we 
need to check whether there exist Mj^. , Mjj- £ Ai{i = 0, . . . ,n) such that Atr.deg Q(Mofco/Moj| 
. . . , Mnk„ /Mnj„) /Q. = n. This can be done with the difference characteristic set method via 
symbolic computation [15]. In this section, a criterion will be given to check whether a 
Laurent difference system is essential in terms of their supports, which is conceptually and 
computationally simpler than the naive approach based on the characteristic set method. 

Let Bi = YYj=iY[k>o(yj''^)'^''^'' ~ 1, . . . ,m) be m Laurent difference monomials. We 
now introduce a new algebraic indeterminate x and let 



Ty, 



sat(R, . . . , Rm), 



supports 



s 
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be univariate polynomials in Z[.t]. If ord(Sj, y^) = — oo, then set dij = 0. The vector {dii,di2, 
. . . , din) is called the symbolic support vector of Bi. The matrix M = {dij)mxn is called the 
symbolic support matrix of Bi, . . . , B„i- 

Note that there is a one-to-one correspondence between Laurent difference monomials 
and their symbolic support vectors, so we will not distinguish these two concepts in case 
there is no confusion. The same is true for a set of Laurent difference monomials and its 
symbolic support matrix. 

Definition 3.16 A matrix M = {dij)mxn over Q{x) is called normal upper-triangular of 
rank r if for each i < r, da and di^i-k = (1 < A; < i — 1), and the last m — r rows are 
zero vectors. 

A normal upper-triangular matrix is of the following form: 

/ ail * • • • * • • • * \ 
022 ■ ■ ■ * ■ ■ ■ * 

• • • Orr • • • * 

•••0 ••• 



\0 •••0 •••0/ 

Definition 3.17 A set of Laurent difference monomials Bi,B2,. . . ,Bjn is said to be in r- 
upper-triangular form if its symbolic support matrix M is a normal upper triangular matrix 
of rank r. 

The following lemma shows that it is easy to compute the difference transcendence degree 
of a set of Laurent difference monomials in upper-triangular form. 

Lemma 3.18 Let Bi, . . . , B^ be a set of Laurent difference monomials in r -upper-triangular 
form. Then Atr.degQ(5i, . . . , B^) /Q = r. 

Proof: Prom the structure of the symbolic support matrix, for i = 1, ... ,r, Bi = nj=i nfc>o(2/j 
with OTd{Bi, yi) > and B^+i = • • • = = L Let B'- = YYj=i Uk>oiyfy'"'- Then 

Atr.degQ(Si,...,B^)/Q 
= Atr.dcgQ(Bi,...,B,)/Q 

> Atr.degQ(yr+i, . . . ,yn){Bi, . . . , Br)/Q{yr+i, ■■■,yn) 
= AtT.degQ{B[,...,B',)/Q. 

So it suffices to prove Atr.dcgQ(5j, . . . , B',) /Q = r. 

If r = 1, i?[ is anonconstant Laurent difference monomial in yi, so Atr.degQ(i?J^)/Q = 1. 
Suppose we have proved for the case r — 1. Let B" = Ylj=iYlk>oiyj'^)'^''^'' ^ then by the 
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hypothesis, Atr.degQ(5i', . . . , B'l 



V-i 



)/Q = r - 1. Thus, 



r > Atv.degq{B[,...,B',)/Q 

= Atr.degQ{S;)/Q + Atr.degQ(S;, . . . , B',)/Q{B',) 

> 1 + AtT.degq{yr){B[, B',_^)/q{yr) 

= 1 + Atr.degQ(S;', . . . , B';_^)/Q = r. 



So Atv.degQ{Bi,...,Bm)/Q 



= r. 



□ 



In the following, we will show that each set of Laurent difference monomials can be 
transformed to an upper-triangular set with the same difference transcendence degree. Here 
we use three types of elementary matrix transformations. For a matrix M over Q[x], Type 
1 operations consist of interchanging two rows of M, say the i-th and j-th rows, denoted by 
r[i, j]; Type 2 operations consist of adding an /(x)-multiple of the j-th row to the i-th row, 
where f{x) £ Q[x], denoted by [i + j{f{x))]; and Type 3 operations consist of interchanging 
two columns, say the i-th and j-th columns, denoted by c[i, j]. In this section, by elementary 
transformations, we mean the above three types of transformations. 

Let Bi, . . . , Bm be Laurent differential monomials and M their symbolic support matrix. 
Then the above three types of elementary transformations of M correspond to certain trans- 
formations of the difference monomials. Indeed, interchanging the i-th and the j-th rows of 
M means interchanging Bi and iJj, and interchanging the i-th and the j-th columns of M 
means interchanging yi and yj in i?i , . . . , B^a (or in the variable order) . Multiplying the i-th 
row of M by a polynomial f{x) = a^x^^ -|- ad-ix'^~^ + • • • + oq € and adding the result 
to the j-th row means changing Bj to 11^=0 • 

Lemma 3.19 Let Bi,...,Bm be Laurent difference monomials and Ci,...,Cm obtained 
by successive elementary transformations defined above. Then Atr.degQ(i?i, . . . , i?m)/Q 



Proof: It suffices to show that Type 2 operations keep the difference transcendence degree. 



That is, for Eto«^^' ^ Q[x], Atr.degQ{B,, B2) /Q = Atr.degQ(5i, nfc=o(^''^i)"'^2)/Q. 

Suppose =Pi/q where pi,q G Z*. Then, clearly, Atr.degQ(i?i)/Q = Atr.degQ(n^^^ 
{a^Bi)P'')/Q. Thus, Atr.degQ(Bi, nLo(^'^i)"'=^2)/Q = Atr.degQ( Ut=oi<^'' BiY' , DLi 



Theorem 3.20 Let Bi, . . . ,Bm be a set of Laurent difference monomials with symbolic sup- 
port matrix M. Then Atr.degQ(5i, . . . , Bm)/Q = rk(M). 

Proof: By Lemma 13.181 and Lemma 13.191 it suffices to show that M can be reduced to a 
normal upper-triangular matrix by performing a series of elementary transformations. 

Suppose M = (dij) Omxn and we denote the new matrix obtained after performing 
elementary transformations also by M. Firstly, perform Type 1 and Type 3 operations when 
necessary to make du 7^ 0. Secondly, try to use du{x) to reduce other elements in the first 
column to by performing Type 2 operations. If there exists an element in the first column 
such that it can not be divisible by du, say d^i) suppose dki{x) = duix)q{x) + r{x) where 
r(x) 7^ and deg(r(x)) < deg(dii(a;)). After performing the transformations [k + l{—q{x))] 



Atr.degQ(Ci, ...,C„)/Q. 




□ 
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and r[l, fc] successively, we obtain a new matrix in which the degree of du strictly decreases. 
Repeat this process when necessary, then after a finite number of steps, we obtain a new 
matrix M such that dii(x) divides each nonzero element in the first column, and by using 
dii (x) to perform Type 2 operations we obtain 

Now we repeat the above process for Mi and whenever Type 3 operations are performed for 
Ml, we assume the same transformations are performed for the whole matrix M. In this 
way, after a finite number of steps, we obtain a normal upper-triangular matrix M. □ 

Example 3.21 Let Bi = yiy2 and B2 = yi^\2'^- Then the symbolic support matrix of Bi 

and B2 is M = ( \ \ ] . Then rk(M) = I I T f Thus, by Theorem [KM. if 
\ X X J I 2 ija ^ b. 

a ^ b, Bi and B2 are transformally independent over Q. Otherwise, they are transformally 

dependent over Q. 

Consider the set of generic Laurent difference polynomials defined in ([T]): 

k 

Pi = UioMiQ + ^ UikMik (i = 1, . . . , m). 

k=l 

Let j3ik be the symbolic support vector of Mik/MiQ. Then the vector Wi = Ylk=o'^'^kf^ik is 
called the symbolic support vector of Pj and the matrix Mp whose rows are wq, . . . ,Wn is 
called the symbolic support matrix of Pq, . . . ,Pn- In terms of Mp, we have the following 
result. 

Theorem 3.22 Follow the above notations. Then the following three conditions are equiv- 
alent. 

1. Pq, . . . ,P„ form a Laurent transformally essential system. 

2. There exist Mik- (i = 0, . . . , n) with 1 < ki < li such that the symbolic support matrix 
of Moko/Moo, . . . , Mnk,jMno is of rank n. 

3. The rank of Mp is equal to n. 

Proof: The equivalence of 1) and 2) is a direct consequence of Theorem 13.201 and Def- 
inition 13.41 And the equivalence of 2) and 3) follows from the fact that det(Mpj) = 

Y.ko,...,U-,k„ nj=ojVj "ifci det((^ofco' • • • , Afc,, • • • , l^nknV), where Mp j is the matrix obtained 
by deleting the {i + l)-th row from M. □ 

We will end this section by introducing a new concept, namely supper-essential systems, 
through which one can identify certain Pj such that their coefficients will not occur in the 
sparse difference resultant. This will lead to the simplification in the computation of the 
resultant. Let T C {0, 1, . . . ,n}. We denote by Pt the Laurent difference polynomial set 
consisting of Pj E T), and Mp^ its symbolic support matrix. For a subset T C {0, 1, . . . , n}, 
if card(T) = rk(Mp^), then Px, or {Ai : i £ T}, is called a transformally independent set. 
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Definition 3.23 Let T C {0, 1, . . . , n}. Then we call T or ¥j supper-essential if the follow- 
ing conditions hold: (1) card(T) — rk(Mp^) = 1 and (2) card(JI) = rk(Mpjj) for each proper 
subset J o/ T. 

Note that supper-essential systems are the difference analogue of essential systems intro- 
duced in [27] and also that of rank essential systems introduced in [23] . Using this definition, 
we have the following property, which is similar to Corollary 1.1 in [27j . 

Theorem 3.24 If {Pq, . . . ,IPn} is a Laurent transformally essential system, then for any 
T C {0, 1, . . . , n}, card(T) — rk(Mp^) < 1 and there exists a unique T which is supper- 
essential. In this case, the sparse difference resultant of Fq, . . . ,Fn involves only the coeffi- 
cients ofFi (ieJ). 

Proof: Since n = rk(Mp) < rk(Mp^)+card(P)-card(PT) = n+l+rk(Mp^)-card(T), we have 
card(T)-rk(Mp^) < 1. Since card(T) -rk(Mp^) > 0, for any T, either card(T) -rk(Mp^) = 
or card(T)— rk(Mp^) = 1. From this fact, it is easy to show the existence of a supper-essential 
set T. For the uniqueness, we assume that there exist two subsets Ti, T2 C {1, . . . , m} which 
are supper-essential. Then, we have 

rk(Mp^^^^J < rk(Mp,J + rk(Mp^J - rk(Mp^^^^J 

= card(Ti) - 1 + card(T2) - 1 - card(Ti n T2) 
= card(Ti UT2) - 2, 

which means that Mp is not of full rank, a contradiction. 

Let T be a supper-essential set. Similar to the proof of Theorem 13. 6^ it is easy to show 
that [Pj]igTnQ{uj}jgTr is of codimension one, which means that the sparse difference resultant 
of Pq, . . . , Pn involves only of coefficients of Pj (i G T). □ 

Using this property, one can determine which polynomial is needed for computing the 
sparse difference resultant, which will eventually reduce the computation complexity. 

Example 3.25 Continue from Example \3.12l It is easy to show that F = {Po,Pi,P2} is 
a Laurent transformally essential system and Po,Pi constitute a supper- essential system. 
Recall that the sparse difference resultant ofF is free from the coefficients 0/P2. 

4 Basic properties of sparse difference resultant 

In this section, we will prove some basic properties for the sparse difference resultant R(uo, 
...,u„). 

4.1 Necessary condition for existence of nonzero solutions 

In this section, we will first give a condition for a system of Laurent difference polynomials to 
have nonzero solutions in terms of sparse difference resultant, and then study the structures 
of nonzero solutions. 

To be more precise, we first introduce some notations. Let • • • ; be a Laurent trans- 
formally essential system of Laurent monomial sets. Each element {Fq, . . . ,Fn) € I^{Ao) x 
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• • • X C{An) can be represented by one and only one point (vg, . . . , v„) G x ■ • • x 
where Vj = {viQ,Vii, . . . ,Vii.) is the coefficient vector of Fi and £ is some difference field 
extension of Q {£ is not fixed but depends on the set Fi). Let Zq^Aq, . . . ,An) be a set con- 
sisting of points (vq, . . . , v„) such that the corresponding Fj = (i = 0, . . . , n) have nonzero 
solutions. That is, 

Zo{Ao, . . . , An) = {(vo, . . . , v„) : Fo = • • • = F„ = have 

a common nonzero solution}. (5) 

The following result shows that the vanishing of sparse differential resultant gives a necessary 
condition for the existence of nonzero solutions. 

Lemma 4.1 Zo{Ao, An) <Z V(Res^o,...,^„) . 

Proof: Let Fq, . . . ,¥n be a generic Laurent transformally essential system corresponding to 
• • • 1 An with coefficient vectors uq, . . . , u„. By Definition l3.81 Res^g^...^^^^ G [Pq, Pi, ... , Pn] 
nQ{uo, . . . ,u„}. For any point (vq, . . . , v„) G Zo(A, ■ ■ ■ let (Pq, • • • ,P„) G_£(A)_x 

• • • X C{An) be the difference polynomial system represented by (vq, . . . , v„). Since Pq, 

have a nonzero common solution, ReSy^o,...,^^ vanishes at (vq, . . . , v„). □ 



Example 4.2 Continue from Example \3.11[ Suppose T = Q(x) and crf{x) = f{x + 1). 
In this example, we have Respp^p^^Pj 7^ 0. But yi = 0, 7/2 = consist of a zero solution of 
Pq = Pi = P2 = 0. This shows that Lemma \4-1\ is not correct if we do not consider nonzero 
solutions. This example also shows why we need to consider nonzero difference solutions, 
or equivalently why we consider Laurent difference polynomials instead of usual difference 
polynomials. 

Remark 4.3 // Problem \3.15\ can be solved positively, then R = also gives a sufficient 
condition for Pq = • • • = P^ = to have a nonzero solution in certain sense. 

The following lemma refiects the structures of the nonzero solutions. 

Lemma 4.4 Let Aq, . . . , An be a Laurent transformally essential system and R = Res^Q^...^yi„ . 
Then there exists some t such that deg(R, Uro) > 0. Suppose ¥i = is a system represented 
by (vo,...,v„) G Zo{Ao,. . . ,An) and ^^(vq, . . . , v„) / 0. If ^ is a common nonzero 
difference solution of Pj = 0(i = 0, . . . , n), then for each k, we have 

Proof: Since X = [N(Po), . . . ,N(P„)] : m is a reflexive prime difference ideal and R G X, 
there exists some r and k such that deg(R, Urfc) > 0. By Lemma [4.8| deg(R, Uto) > 0. 
Denote N(Pi) = MjP^ {i = 0, . . . , n). For each j = 1, . . . , /q, by equation ^ with /c = 0, 
the polynomial ^MrMrj - ^M^M^o e ^- Thus, if ^ is a common nonzero difference 
solution of Pi = 0, then |^(vo, . . . , v„) • Mrj{C) - ^(vq, • • . , v„)M,o(0 = 0. Since 
|^(vo,...,v„) /O, dSD follows. ' □ 
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4.2 Sparse difference resultant is transformally homogeneous 

We now introduce the concept of transformally homogeneous polynomials. 

Definition 4.5 A difference polynomial f G J~{yo, ■ ■ ■ lUn} is called transformally homoge- 
neous if for a new difference indeterminate X, there exists a difference monomial M{X) in X 
such that f{Xyo, . . . , Xyn) = M{X)p{yo, . . . , //deg(M(A)) = m, f is called transformally 
homogeneous of degree m. 

The difference analogue of Euler's theorem related to homogeneous polynomials is valid. 



Lemma 4.6 / G J-{yQ,yi, . . . ,yn} is transformally homogeneous if and only if for each 
r G No, there exists mr G Nq such that 

^ (r) df{yo,...,yn) _ f 

i=Q dy\ 

Proof: Denote Y = {yo, . . . ,yn) temporarily. Suppose / is transformally homoge- 

neous. That is, there exists a difference monomial M(A) = nr=o('^^'''')'"'^ such that /(AY) = 
M(A)/(Y). Then j;;-^, ,M-2^(AY) = -|,(AY)»!|g#2 = = HMmm = 

j^f{XY). Substitute A = 1 into the above equality, we have X]"=o ^i^^^^w ~ "^r/- 

Suppose ord(/,Y) = tq. Then for each r < tq, A^^^ = A^ Etovi'^ ^i>^^) = 
^"^o(Ayi)(')^(AY) = m^/(AY). So /(AY) is homogeneous of degree nir in A^^'). Thus, 

/(AY) = fiXyo, . . . , Xyn, X^'^y^^\ . . . , A(i)y«; . . . ; X^^<^hi['\ . . . , A('^«)y^)) = KU^^'^r^ fin 
Thus, / is transformally homogeneous. □ 

Sparse difference resultants have the following property. 

Theorem 4.7 The sparse difference resultant is transformally homogeneous in each Uj which 
is the coefficient set o/Pj. 

Proof: Suppose ord(R, Uj) = hi > 0. Follow the notations used in Theorem 13. 6i By 
Lemma 13.91 R-(u; Co; • • • ) Cn) = 0. Differentiating this identity w.r.t. u\j^ (j = 1, . . . , k) 
respectively, we have 



+^(_^|iM)"=)=0. (7) 



In the above equations, r;.. {k = 0, . . . ,hi;j = 0, . . . ,li) are obtained by replacing mo by 

9« -J 



Q {i = 0,1, . . . ,n) in each 7^, respectively. 



Multiplying ([7]) by u- • and for j from Itok, adding them together, we get 2^7=1 ^ii (!■) + 
-^Cf^ = 0. Thus, the difference polynomial fk = E'/=o4f vanishes at (Co, • • • ,Cn)- 
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Since ord(/fc,njo) < ord(R, Ujo) and deg(/fc) = deg(R), by Lemma [3.9| there exists an 
mfc € Z such that fk = m^R. Thus, by Lemma |4.6| R is transformahy homogeneous in Uj. 

□ 

4.3 Order bound in terms of Jacobi number 

In this section, we wih give an order bound for the sparse difference resultant in terms of 
the Jacobi number of the given system similar to the differential case. 

Consider a generic Laurent transformahy essential system {Po,...,Pn} defined in ([1]) 
with Uj = {uio,Uii, . . . , lij/J being the coefficient vector of Pj = 0, . . . , n). Suppose R is the 
sparse difference resultant of Pq, . . . ,Pn- Denote ord(R, Uj) to be the maximal order of R 
in Uik (/c = 0, . . . , /j), that is, ord(R, Uj) = maxfeord(R, Uj^). If Uj does not occur in R, then 
set ord(R, Uj) = — oo. Firstly, we have the following result. 

Lemma 4.8 For fixed i and s, if there exists ko such that deg(R, u^^^^) > 0, then for all 

A; G {0, 1, ... , li}, deg(R, u^^jj) > 0. In particular, if ord(R, Uj) = /ij > 0, then ord(R, Uik) = 
hi{k = 0,. . . ,li). 

Proof: Firstly, for each k € {!,..., ^j}, by differentiating R(u; i^Q) • • • i Cn) = w.r.t. u^^\ 
we have -^(u, Co, • • • , Cn) + f^(u, Co, • • • , Cn) ( - feM) = 0. If ko = 0, then is a 

nonzero difference polynomial not vanishing at (u, Co, • • • , Cm) by lemma 13.91 So -^%) ^ 0. 
Thus, deg(R,ttJf ) > for each k. If ko + 0, then ^(u, Co, • • • , Cn) / and ^ / 

follows. So by the case ko = 0, for all k, deg(R,M5^^) > 0. 

In particular, if ord(R, Uj) = /ij > 0, then there exists some ko such that deg(R, u]^^^^ > 
0. Thus, for each /c = 0, . . . , /j, deg(R, ufj^^^) > and ord(R, Uik) = hi follows. □ 

Let A = {oij) be an n X n matrix where Oij is an integer or — oo. A diagonal sum of A 
is any sum aio-(i) + a2(T{2) + • • • + CLna{n) with a a permutation of 1, . . . , n. If A is an m x n 
matrix with M = min{m, n}, then a diagonal sum of j4 is a diagonal sum of any M x M 
submatrix of A. The Jacobi number of a matrix A is the maximal diagonal sum of A, denoted 
by Jac(A). 

Let ord(N(Pj),2/j) = Sij {i = 0, = and ord(N(Pj)) = Sj. We call the 

(n+1) xn matrix A = (sjj) the order matrix oiFo, . . . , P„. By A--, we mean the submatrix of A 
obtained by deleting the (i+l)-th row from A. We use P to denote the set {N(Po), . . . , N(P„)} 
and by P., we mean the set P\{N(Pj)}. We call Jj = Jac(A^) the Jacobi number of the system 
P^, also denoted by Jac(Pj). Before giving an order bound for sparse difference resultant in 
terms of the Jacobi numbers, we first list several lemmas. 

Given a vector 1^ = {ko, ki, . . . , kn) G ^>o^, we can obtain a prolongation of P: 

p[^] = y N(Pj)['='l (8) 

i=0 



18 



Let tj = max{soj +kQ, sij + ki, . . . , Snj + kn}- Then P'-^l is contained in Q[u["^l , Y'-^I] , where 

ut^l = Uf^ouf"' and Yl^l = U]^^y^-'\ 

Denote i/(pl"^l) to be the number of Y and their transforms appearing effectively in 

p["^]. In order to derive a difference relation among Uj (i = 0, . . . , n) from a sufficient 
condition is _^ 

|p[^]| > i/(p[-^]) + l. (9) 

Note that z^(p[^l) < |Y[^1| = Ei=i(*j + !)• Thus, if [p[^l| > YI^] + 1, or equivalently, 

n 

ko + ki-\ h /cn > ^ max(soj + ko, sij + ki,.. . , Snj + kn) (10) 

is satisfied, then so is the inequality 

Lemma 4.9 Let W he a Laurent transformally essential system and it = (/cq, ki, . . . , k^) G 
Z"^"*^ a vector satisfying il(J\). Then ord(R, Uj) < ki for each i = 0, . . . ,n. 

Proof: Denote to be the set of all monomials in variables YI-^1. Let I = (pl^l) : 

be an ideal in the polynomial ring Q[Y[-^], ut-^l]. Denote U = u[-^l\ U"^q u^^\ Let Q = 
-(I]fc=i UikMik/Mio)^''^ fov i = 0,1,. .. ,n;l = 0,1, ... ,ki. Denote C = {U, (oko, • • • , Coo, • • • , 
Crifc„) • • • 1 Cno)- It is easy to show that (Y^-^], is a generic point of Z. Let Xi = Xn Q[u[-^1]. 
Then Xi is a prime ideal with a generic point C. Since Q(C) C Q{Y^^\U), Codim(Xi) = 
|C/| + ^;Lo(^i + l)-ti-.degQ(C)/Q > |C/| + |p[^]|-tr.degQ(Y[^],[/)/Q = |p[^] | - |y[^1 | > 1. 
Thus, Xi 7^ {0}. Suppose / is a nonzero polynomial in Xi. Clearly, ord(/, Uj) < ki and 
/ E [P] : m n Q{uo, . . . , u„}. By Lemma 13.91 and Lemma 14.81 ord(R, Uj) < ord(/, Uj) < k^. 

□ 



Lemma 4.10 J2&, Lemma 5.6] Let ¥ he a system with Ji > for each i = 0,. . . ,n. Then 
ki = Ji{i = 0, . . . ,n) satisfy [W\) in the equality case. 

Corollary 4.11 Let P 6e a Laurent transformally essential system and Ji > for each 
i = 0, . . . ,n. Then ord(R, Uj) < Jj (i = 0, . . . , n). 

Proof: It is a direct consequence of Lemma 14.91 and Lemma 14.101 □ 

The above corollary shows that when all the Jacobi numbers are not less that 0, then 
Jacobi numbers are order bounds for the sparse difference resultant. In the following, we deal 
with the remaining case when some Ji = —oo. To this end, two more lemmas are needed. 

Lemma 4.12 |^ \20f Let A he an mx n matrix whose entries are 's and 1 's. Let Jac(A) = 
J < min{m, n}. Then A contains an a x b zero suh-matrix with a + b = m + n — J. 
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Lemma 4.13 Let¥ be a Laurent transformally essential system with the following (n+1) xn 
order matrix 

f ^11 i-oo)rxt 

where r + t > n + 1. Then r + t = n + 1 and Jac(A22) > 0. Moreover, when regarded 
as difference polynomials in yi, . . . , yr-i, {^0, ■ ■ ■ , Pr-i} is a Laurent transformally essential 
system. 

Proof: The proof is similar to [23^ Lemma 5.9]. □ 

Theorem 4.14 Let ¥ be a Laurent transformally essential system and R the sparse differ- 
ence resultant o/P. Then 



ord(R, Uj) 



—00 if Ji = —00, 

hi < Ji if Ji > 0. 



Proof: Corollary 14.111 proves the case when Ji > for each i. Now suppose there exists at 
least one i such that Jj = —00. Without loss of generality, we assume = — 00 and let 

= {sij)o<i<n-i;i<j<n be the Order matrix of F^. By Lemma r4.12[ we can assume that 
is of the following form 

An {-oo)rxt 



' A21 A22 

where r + t > n + 1. Then the order matrix of P is equal to 

( ^11 {-oo)rxt 
V ^21 A22 

Since P is Laurent transformally essential, by Lemma [4.13l r+t = n+1 and Jac(7422) > 0. 
Moreover, considered as difference polynomials in yi, . . . , yr-i, P = {po, ■ ■ ■ ,Pr-i} is Laurent 
transformally essential and An is its order matrix. Let Ji = Jac((Aii).). By applying the 
above procedure when necessary, we can suppose that Jj > for each i = 0, ...,r — 1. 
Since [P] fl Q{up) • • • > u„} = [P] n Q{uo, . . . , u^-i}, R is also the sparse difference resultant 
of the system P and u^, . . . , u„ will not occur in R. By Corollary 14. IH ord(R, Uj) < Jj. 
Since Ji = Jac(^22) + Jj > Jj for < i < r — 1, ord(R, Uj) < Jj for < i < r — 1 and 
ord(R, Uj) = —00 for i = r, . . . , n. □ 

Example 4.15 Let n = 2 and Pj has the form 

Po = uoo + uoiyiyf^ Pi = ""10 + unyi, P2 = -"10 + uny2 ^- 

/ 1 -oo\ 

In this example, the order matrix of ¥ is A = I —00 I . Thus Jq = 1, Ji = 2, J2 = 

V -00 1 / 

—00. ^ndord(R, uq) = < Jo,ord(R, ui) = 1 < Ji,ord(R, U2) = —00. 

Corollary 4.16 Let P be supper-essential. Then Ji > for i = 0, . . . , n and ord(R, Uj) < Jj . 
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Proof: From the proof of Theorem 14.141 if Jj = — oo for some i, then P contains a proper 
transformally essential subsystem, which contradicts to Theorem 13. 241 Therefore, Jj > for 
i = 0, . . . ,n. □ 

We conclude this section by giving two improved order bounds based on the Jacobi bound 
given in Theorem 14.141 

For each j £ {1, . . . , n}, let Oj = mm{k G No| 3 z s.t. deg(N(Pj), yj^^) > 0}. In other words, 

Oj is the smallest number such that y^"-'^ occurs in {N(Po)i • • • iN(P„)}. Let B = {sij — Oj) 
be an (n + 1) X n matrix. We call Ji = Jac(-Bj) the modified Jacobi number of the system 
P^-. Denote 7 = X]j=i^j- Clearly, Ji = Ji — 7. Then we have the following result. 

Theorem 4.17 Let ¥ be a Laurent transformally essential system and R the sparse differ- 
ence resultant o/P. Then 



ord(R, Uj 



— 00 if Ji = —00, 

hi < Ji — J if Ji > 0. 



Proof: The proof is similar to [231 Theorem 5.13]. □ 

Now, we assume that P is a Laurent transformally essential system which is not supper- 
essential. Let R be the sparse difference resultant of P. We will give a better order bound 
for R. By Theorem 13.241 P contains a unique supper-essential sub-system Fj. Without loss 
of generality, suppose T = {0, . . . , r} with r < n. Let Aj be the order matrix of Pt and for 
i = 0, . . . , r, let (^t)^ be the matrix obtained from Af by deleting the (i + l)-th row. Note 
that (^t)^ is an r X n matrix. Then we have the following result. 



Theorem 4.18 With the above notations, we have 
ord(R, Uj) = 

Proof: Similarly to the proof of ^23^ Theorem 5.16], it can be proved. □ 



hi < Jac((^T)i) i = 0,... ,r, 
—00 i = r + 1, . . . , n. 



Example 4.19 Continue from Example \4-15\ In this example, T = {0,1}. Then Aj = 
Thus Jac{{AT)^) = 0, Jac{{Aj)^) = 1. 74n(i ord(R, ug) = = Jac((A']r)g), ord(R, ui^ 





1 = Jac((^T)i), ord(R, U2) = —00 



5 A single exponential algorithm to compute the sparse dif- 
ference resultant 

In this section, we give an algorithm to compute the sparse difference resultant for a Laurent 
transformally essential system with single exponential complexity. The idea is to estimate 
the degree bounds for the resultant and then to use linear algebra to find the coefficients of 
the resultant. 
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5.1 Degree bound for sparse difference resultant 

In this section, we give an upper bound for the degree of the sparse difference resultant, 
which will be crucial to our algorithm to compute the sparse resultant. Before proposing the 
main theorem, we first give some algebraic results which will be needed in the proof. 

Lemma 5.1 123\ Theorem 6.2] Let X he a prime ideal in K[xi, . . . ,Xri\ and = X n 
K[xi^ . . . , Xk] for any I < k < n. Then deg(Ifc) < deg(I). 

Lemma 5.2 J28[ Corollary 2.28] Let Vi, . . . , 14- C (r > 2) he pure dimensional projective 
varieties in P^. Then 

r 

i=l C 

where C runs through all irreducihle components of Vi D ■ ■ ■ CiVr. 
Now we are ready to give the main theorem of this section. 

Theorem 5.3 Let Pq, . . . ,Pn he a Laurent transformally essential system of form ([ip with 
ord(N(Pj)) = Si ond deg(N(Pj), Y) = rui. Suppose N(Pj) = Y^^l^^UikNik and Ji is the Jacohi 
numher of {N(Po), . . . , N(P„)}\{N(Pi)}. Denote m = maxi{mi}. Let R(uo, . . . ,u„) be the 
sparse difference resultant o/Pi (i = 0, . . . , n). Suppose ord(R, Uj) = hi for each i. Then the 
following assertions hold: 

1) deg(R) < nr=o("^« + 1)'^''^^ < (m + l)EiLo(-^i+i), where m = maxj{mj}. 

2) R has a representation 

n n n hi 

n n {N^^r^^^^ • R = E E G..N(p,)('=) (11) 

i=0 fc=0 i=0 fc=0 

where G,k G Qrf°', . . . , uL''"^, Y^] and h = maxj/ii + ej such that deg(GifeN(Pi)('=)) < 
[m + l + Y.l=oih^ + l)deg(7Vio)]deg(R). 

Proof: In R, let Ujo be replaced by (N(Pj) — UikNik) /NiQ for each i = 0, . . . , n and let R 
be expanded as a difference polynomial in N(Pj) and their transforms. Then there exist G 

N and polynomials dk such that Uto TlfcLo (A^ifV'^R = ELo EfcLo G',fcN(P,)('=) +r with 
T € Q{u, Y} free from mo- Since T el= [N(Po), . . . , N(P„)] : m, T vanishes identically, for 
InQ{u,Y} = {0} by Theorem ESI Thus, 

n hi n hi 

nn(<V"R=EEG^^N(p,)('=). 

i=Ofc=0 j=0 fc=0 

1) Let J = (N(Po)[^«],...,N(P„)['^"]) : ml'^l be an algebraic ideal in 7^ = Q[Y['^1, uf 

...,u|J'"^] where h = maxj{/ij + Sj} and im^^^ is the set of all monomials in Y'^^. Then 
R G i7 by the above equality. Let rj = (rji, . . . ,r]n) be a generic point of [0] over Q(u) 
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and denote Ci = ~ Ylk=i "^ik Nfo(ri) = 0, . . . , n). It is easy to show that J" is a prime ideal 
in 7^ with a generic point (r/t'^l; S, ci''"', • • • , Cn''"^ and J n Q[uf . . . , uL''"'] = (R), where 
u = Uju[''''\{n|Q''}. Let Hik be the homogeneous polynomial corresponding to N(Pj)('^) 
with xo the variable of homogeneity. Then J'^ = {{Hik)i<i<n-o<k<hJ : m is a prime ideal 
, Ug . . . , uli where m is the whole set of monomials in Y^'''] and xq. And 
deg(jO) = deg(J). 

Since Y{{H,k)i<i<n;0<k<hJ = YiJ'^) U Y{Hik,xo)\JUj,iY{H,k,yf), V(jO) is an irre- 
ducible component oiY {{Hif:)i<i<n-fi<k<h^ 

). By LemmaO deg(jO) < ]Xl=oUkloi^^+^) = 
m=oi^^ + Thus, deg(J) < nr=oK + 1)'^'+'- Since J n Qrf "J, . . . , uL'^-^] = (R), 

by LemmaEIIl deg(R) < deg(J) < fULoC™-* + ^ f"'^^ < + l)Er=o(-^»+i) follows. The last 
inequality holds because hi < Ji by Theorem 14.171 

2) To obtain the degree bounds for the above representation of R, that is, to estimate 
deg(GjfcN(Pj)('^)) and Ojfc, we take each monomial M in R and substitute Uio by (N(Pj) — 
Ea;=i '^ik^ik) l^io ii^to M and then expand it. To be more precise, we take one monomial 

M(u;noo,..., uno) = u^nr=ontoK?^)'" ^ith I7I + ELo Eto = deg(R) for an 
example, where u'^ represents a difference monomial in u and their transforms with exponent 
vector 7. Then 

n hi li ^ n hi 

M(u; .00, . . . , u„o) = u'' n n - E ^^>^N,>^f^) "7 H 11 ■ 

i=Ok=a fe=l i=Qk=0 

When expanded, every term of Ui=oUk=o {^io^)'^"' is of degree bounded by I7I + 
Er=o EkLoi^i + '^)dik < (m + l)deg(R) in uf"', . . . , ut'^ and Y^^\ Suppose R = J2m ^mM 
and ttjfc > maxA/jdjfc}. Then 

n hi n hi 

j=0 A:=0 i=0 fc=0 

with deg(G,fcN(Pi)(^-)) < (m + l)deg(R) + ELo EfcLo deg(iVio)aifc. Clearly, we can take 
aik = deg(R) and then deg(GifcN(Pi)(^')) < + 1 + YJi^oih + l)deg(iVio))deg(R). Thus, 
(fTTj) follows. □ 

For a transformally essential difference polynomial system with degree terms, the second 
part of Theorem 15.31 can be improved as follows. 

Corollary 5.4 Let Pj = Ujo + Ea;=i '^ik^ik {i = 0, . . . ,n) be a transformally essential differ- 
ence polynomial system with m = maxj{deg(Pi, Y)} and Ji the Jacobi number o/{Po, . . . ,Pn}\{Pi}- 
Let R(uo, . . . , u„) be the sparse difference resultant o/Pj (z = 0, . . . , n). Suppose ord(R, Uj) = 
hi for each i and h = max{/ij + Si}. Then R has a representation 

n hi 

R(uo,...,u„) = EE^^i^?^ 

i=0 j=0 

where Gij G Q[uf°^, . . . , u[f"l, YI'*]] such that deg{Gij¥[^^) < (m+l)deg(R) < (m+l)E"=oW+i)+i. 
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Proof: It is direct consequence of Theorem 15.31 by setting A'jo = 1- D 
The fohowing result gives an effective difference Nullstehensatz under certain conditions. 



Corollary 5.5 Let /o, . . . , S J'{yi, ■ ■ ■ , Un} have no common solutions with deg(/j) < m. 
Let Jac({/o, . . . , fn}\{fi}) = Ji- If the sparse difference resultant of fo, . . . , fn is nonzero, 
then there exist Hij G T{yi, ... ,yn} s.t. Y17=o Y^f=o ^ij fi^'' = 1 "-''^d deg{Hijf^^^) < (m + 

i)Er=o(-^»+i)+i. 

Proof: The hypothesis imphes that P(/i) form a transformahy essential system. Clearly, 
R(uo, . . . ,Un) has the property stated in Corollary 15.41 where Uj are coefficients of IP(/i). 
The result follows directly from Corollary 15.41 by specializing Uj to the coefficients of /j. □ 

5.2 A single exponential algorithm to compute sparse difference resultant 

If a polynomial R is the linear combination of some known polynomials Fi{i = 1, . . . , s), that 
is i? = LLiFi, and we know the upper bounds of the degrees of R and HiFi, then a 

general idea to estimate the computational complexity of R is to use linear algebra to find 
the coefficients of R. 

For sparse difference resultant, we already have given its degree bound and the degrees 
of the expressions in the linear combination in Theorem 15.31 

Now, we give the algorithm SDResultant to compute sparse difference resultants based 
on the linear algebra techniques. The algorithm works adaptively by searching for R with 
an order vector (/iq, . . . , hn) € Nq"*"^ with hi < Ji by Theorem 15.31 Denote a = Y17=o 
We start with o = 0. And for this o, choose one vector (ho, . . . ,hn) at a time. For this 
{ho, . . . ,hn), we search for R from degree d = 1. If we cannot find an R with such a degree, 
then we repeat the procedure with degree d-\-l until d > nr=o('^* l)'^'^^- that case, we 
choose another {ho, . . . , with Y17=o ^« ~ -^^^ (^Oi • • • j hn) with hi < Ji and 

Y17=o hi = o, a cannot be found, then we repeat the procedure with o + 1. In this way, we 
will find an R with the smallest order satisfying equation (llip . which is the sparse resultant. 

Theorem 5.6 Let ¥q, . . . ,P„ be a Laurent transformally essential system of form (Cp. De- 
note P = {N(Po), . . . ,N(P„)}, Ji = Jac(Pj), J = maxjJi and m = max^^odeg(Pi, Y). Algo- 
rithm SDResultant computes sparse difference resultant R o/Pq, . . . ,Pn with the following 
complexities: 

1) In terms of the degree hound D ofH, the algorithm needs at most 0{D^^^'^\nJ)^^'"^^) 
'^-arithmetic operations, where I = 'Yl!i=oi^i + 1) ■^^■^^ "^i- 

2) The algorithm needs at most 0(m'^^"'"'^^(nJ)*~^^''^)) Q-arithmetic operations. 

Proof: The algorithm finds a difference polynomial P in Q{uo, . . . ,u„} satisfying equation 
(llip . which has the smallest order and the smallest degree in those with the same order. 
Existence for such a difference polynomial is guaranteed by Theorem 15.31 By the definition 
of sparse difference resultant, P must be R. 

We will estimate the complexity of the algorithm below. Denote D to be the degree 
bound of R. By Theorem 15.31 D < (m + l)J^i=o('^!+i). In each loop of Step 3, the complexity 
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Algorithm 1 — SDResultant(Po, ■ ■ ■ 



Input: A generic Laurent transformally essential system Pq, . . . ,P„,. 
Output: The sparse difference resultant R(uo, . . . , u„) of Pq, . . . , Pn- 

1. For i = 0, . . . ,n, set N(Pi) = Ek=o^ikNik with deg(iVio) < deg(Aife). 
Set rrii = deg(N(Pi)), rriio = deg(Ajo), Uj = coeff(Pi) and |uj| = + 
Set Sij = ord(N(Pj), yj), A = (sij) and compute Jj = Jac{A:^. 

2. Set R = 0, o = 0, m = maxj{mj}. 

3. While R = do 

3.1. For each {ho, . . . , hn) G Nq"*'"'^ with X^ILo ^ — "^i 

3.1.1. U = U^=ouf ^ = maxi{/ii + a}, d = I. 

3.1.2. While R = and d < ]Xl=oi^i + 1)^"^^ do 

3.1.2.1. Set Ro to be a homogeneous GPol of degree d in U. 

3.1.2.2. Set Co = coeff(Ro, U). 

3.1.2.3. Set Hij{i = 0, . . . ,n; j = 0, . . . , hi) to be GPols of degree 
[m + l + Er=o(/i* + ^)^io]d - - 1 in Y^^^lU. 

3.1.2.4. Set = coeff(//ij, U U). 

3.1.2.5. Set r to be the set of coefficients of IlILo TlfcLol^io VR-Q- 
EILo E}=o ^*j(N(P»))(^) as a polynomial in Y^^\U. 

3.1.2.6. Solve the linear equation P = in variables cq and Cij. 

3.1.2.7. If Co has a nonzero solution, then substitute it into Ro to 
get R and go to Step 4, else R = 0. 

3.1.2.8. d:=d+l. 

3.2. o:=o+l. 

4. Return R. 

/*/ GPol stands for generic algebraic polynomial. 

/*/ coeff(P, V) returns the set of coefficients of P as an ordinary polynomial in variables V. 



of the algorithm is clearly dominated by Step 3.1.2, where we need to solve a system of linear 
equations V = over Q in co and Cjj. It is easy to show that |co| = C^^^^^) ~ 
(''"Xnl^i?''^'^ ^l^ere L = Y:^=oih^ + l){k + 1) ^ud di = [m + I + Y:^=oih^ + l)nl^o\d. 
Then P = is a linear equation system with N = C^l^Y^) + T^i=o{h^ + 1) {^''""l^^^^ti)^^^^) 
variables and M = ('^^2+r!(fc+i)^^) equations. To solve it, we need at most (max{M, A^})"^ 
arithmetic operations over Q, where uj is the matrix multiplication exponent and the cur- 
rently best known uj is 2.376. 

The iteration in Step 3.1.2 may go through 1 to nr=o("^« + 1)'^''^^ < (m + l)EILo(-^«+i), 
and the iteration in Step 3.1 at most will repeat nr=o("^» + ^) — {n + 1){J + 1) times, where 
J = maxjjj. And by Theorem 15.31 Step 3 may loop from o = to Yll=o{Ji + The whole 
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algorithm needs at most 

0=0 hi<J^ d=l 

arithmetic operations over Q. In the above inequahties, we assume that (m+l)S«=o(-^»+i)+i > 
l{n + 1)J and / > (n + 1)^, where I = Y^^^^ih + !)• Our complexity assumes an 0(1)- 
complexity cost for all field operations over Q. Thus, the complexity follows. □ 



Remark 5.7 As we indicated at the end of Section 3.3, if we first compute the suuper- 
essential set T, then the algorithm can be improved by only considering the Laurent difference 
polynomials Pj (i G T) in the linear combination of the sparse resultant. 

Remark 5.8 Algorithm SDResultant can be improved by using a better search strategy. If 
d is not big enough, instead of checking d+1, we can check 2d. Repeating this procedure, we 
may find a k such that 2^ < deg(R) < 2^^^. We then bisecting the interval [2'^, 2'^"'"-^] again 
to find the proper degree for R. This will lead to a better complexity, which is still single 
exponential. 



For difference polynomials with non-vanishing degree terms, a better degree bound is 
given in Corollarv 15.41 Based on this bound, we can simplify the Algorithm SDResultant 
to compute the sparse difference resultant by removing the computation for N(Pj) and NiQ 
in the first step where N^q is exactly equal to 1. 

Theorem 5.9 Algorithm SDResultant computes sparse difference resultants for a trans- 
formally essential system of the formFi = Uio+Y^^'/l^iUikNik with at most 0{n^'^'^^ J^^^'^mP^^'"^ 
'^-arithmetic operations. 

Proof: Follow the proof process of Theorem 15.61 it can be shown that the complexity is 

(9(-^3.376 jO(n)^0(n/j2))_ □ 

6 Difference resultant 

In this section, we introduce the notion of difference resultant and prove its basic properties. 

Definition 6.1 Let mis^r be the set of all difference monomials in Y of order < s and degree 
< r. Let u = {uMjMemsr be a set of difference indeterminates over Q. Then, 

is called a generic difference polynomial of order s and degree r. 
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Throughout this section, a generic difference polynomial is assumed to be of degree 
greater than zero. Let 



Pi = Uio + 



E 



(¥Wr(i = o,i...., 



n) 



(12) 




1 < |q:| < rrii 



be generic Laurent difference polynomials of order Sj, degree mi, and coefficients Uj respec- 
tively. Since {1, yi, . . . , y„} is contained in the support of each Pj, clearly, they form a 
supper-essential system and the sparse difference resultant Resp(,^...^p^(uo, . . . , u„) exists. We 
define RespQ^...^p„(uo, . . . , u„) to be the difference resultant of Pq, . . . ,Pn. That is. 

Definition 6.2 Let Pj (i = 0, 1, . . . , n) be a generic difference polynomial system of the form 
Then the difference resultant R(uo, . . . , u„) o/Pq, . . . ,P„ is defined as the irreducible 
difference polynomial contained in [Pq, . . . , Pn] H Q{uo, . . . , u„} of minimal order in each Uj, 
which is unique up to a factor in Q. 

Difference resultants hold all the properties we have proved for sparse difference resultants 
in previous sections. Apart from these, in the following, we will show difference resultants 
possess other better properties. Firstly, we will give the precise degree for the difference 
resultant, which is of BKK-type [HE]. Before doing so, we need some results about algebraic 
sparse resultants. 

Let )C[X] = K[ xi, . . . , be the polynomial ring defined over a field K,. For any vector 
a = (ai, . . . , a„) G Z", denote the Laurent monomial x\^X2^ ■ ■ ■ by X". Let Bq, ■ ■ ■ ,13n C 
Z" be subsets which jointly span the affine lattice Z". Suppose = {0, . . . ,0) £ Bi for each 
i and \Bi\ = li + 1 > 2. Let 



be generic sparse Laurent polynomials defined w.r.t Bi {i = 0,1, ... , n). Bi or {X" : a € Bi] 
are called the support of Fj. Denote Cj = {cia)a&Bi ^ — Ui(cj\{cio}). Let Qi be the 
convex hull of Bi in M", which is the smallest convex set containing Bi. Qi is also called the 
Newton polytope of Fj, denoted by NP(Fj). In [27], Sturmfels gave the definition of algebraic 
essential set and proved that a necessary and sufficient condition for the existence of sparse 
resultant is that there exists a unique subset i which is essential. Now, we restate the 

definition of essential sets in our words for the sake of later use. 

Definition 6.3 Follow the notations introduced above. 

• A collection of {Bi}i^ j, or {Fjjjgj of the form IW\) . is said to be algebraically inde- 
pendent i/ tr.degQ(c)(Fj — Cio : i G J)/Q(c) = |J|. Otherwise, they are said to be 
algebraically dependent. 

• A collection of {Bi}i,z / is said to be essential if {Bi}i,:j is algebraically dependent and 
for each proper subset J of I, {Bi}i<z j are algebraically independent. 





(13) 



aeeA{o} 
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In the case that {BQ,...,Bn} is essential, the degree of the sparse resultant can be 
described by mixed volumes. 

Theorem 6.4 (^7]) Suppose that {Bq, . . . , Bn} is essential. For each i E {0, 1, . . . , n}, the 
degree of the sparse resultant in Cj is a positive integer, equal to the mixed volume 

MiQo,...,Q^-l,Qi+l,■■■,Qn)= J2 (-ir-l-'lvoK^ Q,) 

JC{0,...,i-l,i+l,...,n} jGJ 

where vol(Q) means the n- dimensional volume of Q d and Qi + Q2 means the Minkowski 
sum of Qi and Q2. 

Now, we give the first main result of this section. 

Theorem 6.5 Let Pj (i = 0, . . . ,n) be generic difference polynomials in Y = {yi, . . . ,yn} 
with order Si, degree rui, and coefficients Uj respectively. Let R(uo, . . . ,\in) be the difference 
resultant o/Pq, . . . , IPn- Denote s = Y17=o ^i- ^^^^^ R-(uo, . . . , u„) is also the algebraic sparse 
resultant 0/ P|j* p[f as polynomials in YW, and for each i G {0,1,..., n} and 

/t 5 ... 5 5 (§2 ; 

ord(R, Uj) = s — Si (14) 
deg(R, uf ^) = M{{Qji)j^ifi<i<s-sj, Qio, Qi,k~i, Qi,k+i, • • • , Qi,s-s,) (15) 

where Qji is the Newton polytope o/P^^ as a polynomial in Y^^^ and u-'^^ = {uf^ , Uia € Uj}. 

Proof: Reg ard Pf ^ {i = 0, . . . ,n; k = 0, . . . , s — Si) as polynomials in the n{s + 1) variables 
"^'^^ = {2/1, • • • , yn, • ■■,yn\---, ui^K ■ • • , Un^}, and we denote its support by Bik- Since 
the coefficients of F^^^ can be treated as algebraic indeterminates, P^'^^ are generic sparse 
polynomials with supports Bn^^ respectively. Now we claim that: 

CI) B = {Bik ■ < i < n;0 < k < s — Si} is an essential set. 

C2) B = {Bik :0<i<n;0<k<s-Si} jointly spans the affine lattice Z"(^+i). 

Note that \B\ = n{s + 1) + 1. To prove CI), it suffices to show that any n{s + 1) of 
distinct P) are algebraically independent. Without loss of generality, we prove that for a 
fixed / G {0, . . . , s — sq}, 

Q — rmC^OA m> p{'-i) raC'+l) Tp(s-so)x 

— ll^i jl<j<n;0<fc<s-Si;Jr'05 • • • 5^0 '^0 i • • • ; ^0 J 

is an algebraically independent set. Clearly, {y^P , • • • , y^^^'^]] = 1, . . . , n} is a subset of the 

(k) (k) / (k)\ 

support of . Now we choose a monomial from each P^ and denote it by m{¥;- ). Let 

mW^ / ^ 0<A;</-1 (fc). j yf^""^ 0<k<so 

m(PX = < I i,^ and m(P ) = < (.-,+k) 
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For each i G {2, . . . , n}, let 



So m{Si) is equal to {yj*^ : 1 < i < ?^}, which are algebraically independent over Q. 
Thus, the n(s + 1) members of 5*; are algebraically independent over Q. For if not, all the 
Pp-* — u-^ G Si) are algebraically dependent over Q(v) where v = U-Lgu|''~'''\{n|Q~ 
Now specialize the coefficient of m(P-'^^) in P^^^* to 1, and all the other coefficients of P^'^'' —u'^ 
to 0, by the algebraic version of Lemma l2.2l {m(p|'^^) : P^^^^ G Si} are algebraically dependent 
over Q, which is a contradiction. Thus, claim CI) is proved. Claim C2) follows from the fact 
that 1 and yI^I are contained in the support of Pg ''^ . 

By CI) and C2), the sparse resultant of (P.p^)o<i<n.;0<fc<s-Si exists and we denote it 
by G. Then (G) = ((Pf ^)o<i<n;0<fc<.-.,) H Q[u!,'"'°l, • • • ,uk^-'"^], and by Theorem El 



deg(G, uf ^) = M{{Qji)j^i^Q<i<s--sp Qio, • • • , Qi,k-i, Qi,k+i, Qi,s-s,), where uf ^ = {u[q\ 

(k) 

. . . , ulJ , . . .). The theorem will be proved if we can show that G = c • R for some c G Q. 
Since G G [Pq, . . . , Pn] and ord(G, Uj) = s — Sj, by Lemma [3^ ord(R, Uj) < s — Si for each 

i = 0, . . . ,n. If for some i, ord(R,Uj) = hi < s-Sj, then R G {{'!Pf^)jjLi-o<k<s-sj,'^i, ■ ■ ■ 

a contradiction to CI). Thus, ord(R, Uj) = s — Si and R G (G). Since R is irreducible, there 

exists some c G Q such that G = c • R. So R is equal to the algebraic sparse resultant of 

p[f"^°',...,pL^"''"'. □ 

As a direct consequence of the above theorem and the determinant representation for 
algebraic sparse resultant given by D' Andrea [8], we have the following result. 

Corollary 6.6 The difference resultant for generic difference polynomials Pj,i = 0, ...,n 
can he written as the form det(Mi)/ det(Mo) where Mi and Mq are matrixes whose elements 
are coefficients of Pj and their transforms up to the order s — Si and Mq is a minor of Mi . 

Based on the matrix representation given in the above corollary, the single exponential 
algorithms given by Canny, Emiris, and Pan [T2] an be used to compute the difference 
resultant. 

Remark 6.7 From the proof of Theorem \6.5\ we can see that for each i and < k < s — st, 
deg(R, u-'^^) > 0. Furthermore, by Lemma [JT^ deg(R, «|q^) > and deg(R, u^^'') > for 
each a. In particular, deg(R, Ujo) > and deg(R, Ujo) > 0. 

Now, we proceed to give a Poisson-type product formula for difference resultant. Let 
li = UjLgUj \ {uoo} and Q(u) be the transformally transcendental extension of Q in the usual 
sense. Let Qg = Q(u)(uoo, • • • , Wqq ^° ^^). Here, Qo is not necessarily a difference overfield of 
Q, for the transforms of uoo are not defined. In the following, we will follow Cohn ^ to obtain 
algebraic extensions Qi of Qo and define transforming operators to make Qi difference fields. 
Consider R as an irreducible algebraic polynomial r(tiQQ ^°^) in Qo[^oo ^°^]- ^ suitable 
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algebraic extension field of Qo) ^(^^oo *°^) ~ ^ *o = deg(r, Uqq '^°^) = deg(R, Uqq *"^) roots 
7i,--- >7to- Thus 

to 

R(uo, ui, ...,Un) = All (4^0-^"^ - 7^ (16) 

T=l 

where A G Qq. Let Zu = Po, • • • , Pn] n Q{uo, . . . , u„}. Then by Definition [621 is an 
essential reflexive prime difference ideal in the decomposition of {R} which is not held by 
any difference polynomial of order less than s — sq in uqq. Suppose R, Ri, R2, ... is a basic 
sequenc^ of R corresponding to lu- That is, lu = U/c>o asat(R, Ri, . . . , Rfc). Regard all 
the Rj as algebraic polynomials over the coefficient field Q(u). Denote ^rO = 7t- Clearly, 
Uqq = 7t-o is a generic point of asat(R). Suppose {i < k) are found in some al- 
gebraic extension field of Qo such that Uqq = 7^^ (0 < i < A;) is a generic point of 
asat(R, Ri, . . . , Rfc). Then let 7r,A;+i be an element such that u^q ''"^'^ = (0 < i < k + 1) 
is a generic point of asat(R, Ri, . . . , R^, Rfc+i). Clearly, 7r,A;+i is also algebraic over Qo- Let 

Ot = Q(u)(woO) • • • ) ftoo 7r J 7t1, • • •)• Clearly, is an algebraic extension of Qo and Qr 

is algebraically isomorphic to the quotient field of Q{uo, . . . , u„}/Xu. Since the quotient field 
of Q{uo, . . . , u„}/Iu is also a difference field, we can introduce a transforming operator CTt- 
into Qr to make it a difference field such that the above isomorphism becomes a difference 
one. That is, cJt-Iqo = o"Iqo and 



cTriuoo) 



u^oo 0<k<s — sq — 1 

lT,k-s-so k> s - So 



In this way, ((7^,0",-) is a difference field. 

Let F be a difference polynomial in Q{uo, ui, . . . , u„} = Q{u, uoo}- For convenience, 

by the symbol F| (s-so)_ , we mean substituting tilf, ""'^'^^ by (t^Jt = 7Tfc {k > 0) into F. 

"00 —"/-^ 

Similarly, by saying F vanishes at Uqq = 7,-, we mean F\ (s-so)_ = 0. The following 

""oo "T'T 

lemma is a direct consequence of the above discussion. 
Lemma 6.8 F G lu if CL^d only if F vanishes at Uqq = 7^. 

Proof: Since lu = Ua;>o asat(R, Ri, . . . , R^) and Uqq = 7ri (0 < i < A;) is a generic 

point of asat(R, Ri, . . . , Hk), the lemma follows. □ 



Remark 6.9 In order to make Qr o difference field, we need to introduce a transform- 
ing operator Gr which is closely related to jr- Since even for a fixed r, generic points of 
asat(R, Ri, . . . , R^) beginning from Uqq ^"'^ = 7,- may not be unique, the definition of Ur also 
may not be unique, which is different from the differential case. In fact, it is a common 
phenomena in difference algebra. Here, we just choose one, for they do not influence the 
following discussions. 

^For the rigorous definition of basic sequence, please refer to [3]. Here, we list its basic properties: i) For 
each k > 0, ord(Rfc , uoo ) = s — sq + k and R, Ri , . . . , Rfc is an irreducible algebraic ascending chain, and ii) 
|Jj.>Q asat(R, Ri, . . . , Rfc) is a reflexive prime difference ideal. 
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Now we give the following Poisson type formula for the difference resultant. 

Theorem 6.10 Let R(uo, . . . , u„) be the difference resultant o/Fq, . . . ,Fn- 

= to- Then there exist S^rk {t = 1, . . . ,tQ; k = 1, . . . ,n) in overfields {Qr-, (Jr) of (Q(u), a) such 

that 

to 

R = A\{Fo{U,---,irn)^'-''\ (17) 

T=l 

so]\ ^,(«-«o)i 



where A G Q(ui, . . . , u„) [ug " x^^oo ° ] ■ ^ote that (11) is formal and should be understood 

in the following precise meaning: Po(^tI) • • • > ^m)^^"^"'^ = *°^ioo+fT~'^''(Sa Uoa{^,T ^°^)°'), 
where (,r = {Cri, ■ ■ ■ lim)- 

Proof: By Theorem 14.71 there exists m G N such that 

5R ^ dn 



E(/XV 
uoa-K = mR. 



00 ^ OUoa 



Setting Uqq = 7t- in both sides of the above equation, we have 



du. 



00 



Oa 



, , =0. 

{s — sn) 
"00 =7t 



Let = (S/£)L(»--o)=^^- Then uoo = -Ea^OaC™ with uJ^q = 7^. That is, 

7, = -fTr^"(Ea^Oae™) = -(E„noae™)(^-^«). Thus, 

to 

R = ^ (noo + ^ uoa6 



r=l 



Suppose Po = ^^00 + Ej=i'"Ojoyj + To. Let = (a!5^/£)|„(.-«o)=^^ 0' = 1, • • • and 

?T = (^ri, • • • ) '^rn)- It remains to show that ^t-q = (^t °')". 

Let Ci = - I]q, Ui« (i = 0, . . . , n). Clearly, C = (u, Co, • • • , Cn) is a generic point of 

= [Po,...,Pn]nQ{uo,..., u„}, where u = UtiuA{n.o}. For each (yM)° = 
byequation©,(YM).=^/^= ft fl ((^^/S^'T"' ^^^^S = £ 



So S ft ft - £ ft ft (W)^^V" ^ e 



'^jfc ^ 

G T.,.. Rv T,emma IfiTRl 

j=l k=0 ^ ^ •'" j=i k=0 ^ ^ 

ft ft i^r-Y'" = (^r Thus, m follows. □ 

j=l fc=0 

Theorem 6.11 The points = {^t1, ■ ■ ■ ^im) {t' = 1, . . . ,to) (fi7| j are generic points of 
the difference ideal [Pi, . . . ,P„] C Q(ui, . . . , u„){Y}. 
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Proof: Clearly, are n-tuples over Q(ui, . . . , u„). For each i = l,...,n, rewrite Pj = 
u^o + Euia n n iyfr^'- Smce C^ = -Euia U n (yf )"^^ and y, = ^J^^, Q + 

a j=l k=l a j=l k=l 

n Si n Si 

E^.« n n ii^JS-Z'T'" = 0. Let a,. = max„a,,. Then n,o U U i&'T' + 

a j=l k=l j=l k=l 

n Si 

E^^aUU ii^J'T'i&'V""' e lu. Thus, by Lemma EH ^r) = u.o + 

a j=l k=l 

E^- n fi(el?r^'' = 0(^ = 1,. ..,n). 

a j=l k=l 

On the other hand, suppose F € Q(ui, . . . ,u„){Y} vanishes at .^t-. Without loss of gen- 
erality, suppose F G Q{ui, . . . , u„, Y}. Clearly, Pi,...,P„ constitute an ascending chain 
in Q{ui, . . . ,u„,Y} with Uiq as leaders. Let G be the difference remainder of F with re- 
spect to this ascending chain. Then G is free from UiQ and F = Gmod [Pi, . . . ,Pn]. Then 
G(^r) = G(u;^T-i, • • • ,&n) = 0, where u = Uf^iUi\{uio} ■ So there exist G N such that 
Gi = Uk ((S)^''T'G(ii;^) G Xu. Thus, Gi vanishes at = Q {i = I, . . . ,n) while |^ 
does not. It follows that G(u;Y) = and F G [Pi,...,Pn]. So ^t- are generic points of 
[Pi,...,P„] cQ(ui,...,u„){Y}. □ 

By Theorem 16. 101 and 16. m we can see that difference resultants have Poisson-type prod- 
uct formula which is similar to their algebraic and differential analogues. 



7 Conclusion and problem 

In this paper, we first introduce the concepts of Laurent difference polynomials and Lau- 
rent transformally essential systems and give a criterion for Laurent transformally essential 
systems in terms of their supports. Then the sparse difference resultant for a Laurent trans- 
formally essential system is defined and its basic properties are proved. Furthermore, order 
and degree bounds for the sparse difference resultant are given. Based on these bounds, an 
algorithm to compute the sparse difference resultant is proposed, which is single exponential 
in terms of the order, the number of variables, and the size of the Laurent transformally es- 
sential system. Besides these, the difference resultant is introduced and its basic properties 
are given, such as its precise order, degree, determinant representation and the Poisson-type 
product formula. 

In the rest of this section, we propose several questions for further study apart form 
Problem 13.151 

It is useful to represent the sparse difference resultant as the quotient of two determinants, 
as done in [U [TT] in the algebraic case. In the difference case. Theorem 16.51 shows that 
difference resultant has such a matrix formula, but for sparse difference resultant, we do not 
have such a formula yet. From (jlip . a natural idea to find a matrix representation is trying 

to define the sparse difference resultant as the algebraic sparse resultant of P = {pf ^ {i = 
0, . . . , n, = 0, . . . , hi)} considered as Laurent polynomials in 

The degree of the algebraic sparse resultant is equal to the mixed volume of certain 
polytopes generated by the supports of the polynomials [23] or [161 P-255]. A similar degree 
bound is given [23^ Theorem 1.3] for the differential resultant. And Theorem 16.51 shows that 
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the degree of difference resultants is exactly of such BKK-type. We conjecture that sparse 
difference resultant has such degree bounds. 

There exist very efficient algorithms to compute algebraic sparse resultants [TOl HH EZl IE] , 
which are based on matrix representations for the resultant. How to apply the principles 
behind these algorithms to compute sparse difference resultants is an important problem. 
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